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ABSTRACT 

RS Ophiuchi began its latest outburst on 2006 February 12. Previous out- 
bursts have indicated that high velocity ejecta interact with a pre-existing red 
giant wind, setting up shock systems analogous to those seen in Supernova Rem- 
nants. However, in the previous outburst in 1985, X-ray observations did not 
commence until 55 days after the initial explosion. Here we report on Swift ob- 
servations covering the first month of the 2006 outburst with the Burst Alert 
(BAT) and X-ray Telescope (XRT) instruments. RS Oph was clearly detected in 
the BAT 14-25 keV band from t = to t ~ 6 days. XRT observationsfrom 0.3- 
10 keV, started at 3.17 days after outburst. The rapidly evolving XRT spectra 
clearly show the presence of both line and continuum emission which can be fit- 
ted by thermal emission from hot gas whose characteristic temperature, overlying 
absorbing column, [A^-jp^/, and resulting unabsorbed total flux decline monoton- 
ically after the first few days. Derived shock velocities are in good agreement 
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with those found from observations at other wavelengths. Similarly, [A'"j|/]vi/ is in 
accord with that expected from the red giant wind ahead of the forward shock. 
We confirm the basic models of the 1985 outburst and conclude that standard 
Phase I remnant evolution terminated by t ~ 10 days and the remnant then 
rapidly evolved to display behaviour characteristic of Phase III. Around t = 26 
days however, a new, luminous and highly variable soft X-ray source began to 
appear whose origin will be explored in a subsequent paper. 

Subject headings: binaries: close - novae, cataclysmic variables - stars: individual 
(RS Oph) stars: supernovae - stars: symbiotic - white dwarfs 



1. Introduction 

RS Ophiuchi is a Symbiotic Recurrent Nova (RN) which had previously undergone 
recorded outbursts in 1898, 1933, 1958, 1967 and 1985 (see Rosino 1987; Rosino & lijima 
1987), with a possible additional outburst in 1907 (Schaeffer 2004). On 2006 February 
12.83UT it was observed to be undergoing a further eruption (Hirosawa 2006), reaching 
magnitude V=4.5 at this time. For the purposes of this paper, we define this as t = to. 
The optical light curve then continued a rapid decline, consistent with that seen in previous 
outbursts (Rosino (1987), AAVSO^). 

The RS Oph binary system comprises a red giant star in a 455.72 ± 0.83 day orbit with 
a white dwarf (WD) of mass near the Chandrasekhar limit (see Dobrzycka & Kenyon 1994; 
Shore et al. 1996; Fekel et al. 2000). Accretion of hydrogen-rich material from the red giant 
onto the WD surface leads to the conditions for a thermonuclear runaway (TNR) in a similar 
fashion to that for Classical Novae (CNe). The much shorter inter-outburst period for this 
type of RN compared to CNe is thought to be due to a combination of the high WD mass 
and a supposed high accretion rate (e.g. Starrfield et al. 1985; Yaron et al. 2005). Such 
models lead to the ejection of somewhat lower masses at higher velocities than those for CN 
models (typically 5000 km s~^ and 10^^ — IO^^Mq respectively for RNe). Spectroscopy of 
RS Oph has indeed shown Ha line emission with FWHM = 3930 km s~^ and FWZl = 7540 
km s~^ on 2006 February 14.2 (t = 1.37 days. Bull (2006)). Superimposed on the broad line 
is an intense and narrow double-peaked structure. 

Unlike CNe, where the mass donor is a low-mass main-sequence star, the presence of 
the red giant in the RS Oph system means that the high velocity ejecta run into a dense 
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circumstellar medium in the form of the red giant wind, setting up a shock system with gas 
temperatures ~ 2.2 x lO^K for vs = 4000 km s~^, where vs is the velocity of the forward 
shock running into the pre-existing wind (see below). Evidence for the presence of such high 
temperature material in outbursts prior to 1985 came from observations of coronal lines in 
optical spectra (Rosino 1987) and the expected deceleration was evidenced by the narrowing 
of the initially broad emission lines (Snijders 1987; Shore et al. 1996). The superimposed 
narrow lines are then from emission and absorption in the red giant wind ahead of the 
forward shock. 

The interstellar absorbing column, Nh = (2.4 ±0.6) x 10^^ cm~^, was determined from 
HI 21cm measurements (Hjellming et al. 1986) and is consistent with the visual extinction 
{E{B — V) = 0.7 ± 0.1) determined from lUE observations in 1985 (Snijders 1987). These 
quantities, together with the observed versus theoretical bolometric luminosity and apparent 
brightness of the red giant, are all consistent with a distance to RS Oph of 1.6 ± 0.3 kpc 
(Bode 1987). 

EXOSAT observations in 1985 (Mason et al. 1987) from t = 55 — 251 days post-outburst 
showed that RS Oph was initially an intense, then a rapidly declining soft X-ray source. Bode 
& Kahn (1985) formulated an analytical model of the evolution of the outburst based on the 
X-ray emission observed by EXOSAT at 55 days, the rapidly increasing radio emission from 
t = 18 days (Padin et al. 1985) and parameters of the red giant wind derived from optical 
observations. They concluded that RS Oph evolved like a Supernova Remnant (SNR), but 
on timescales around 10^ times faster (see below). Subsequently, O'Brien & Kahn (1987) and 
O'Brien et al. (1992) constructed detailed analytical and numerical models of the interaction 
of the ejecta with the circumstellar medium which led to consistent estimates of the outburst 
energy and ejected mass of 1.1 x 10'^^ erg and 1.1 x 10^^ Mq respectively, with the ratio of 
mass loss rate to outflow speed in the red giant wind being estimated as 6 x 10^^ g cm"^, 
compatible with that from isolated red giants. Perhaps not surprisingly, this relatively simple 
model failed to agree fully with the observed spectral evolution of the X-ray emission. In 
particular, it was difficult to reconcile both the low-energy and higher energy behaviour of 
the emission (see also Itoh & Hachisu (1990) who primarily explored the interaction of shells 
from two separate outbursts). 



2. Swift Observations 

Swift (Gehrels et al. 2004) is a multiwavelength mission primarily designed to detect 
gamma-ray bursts (GRBs). However the wide- field, hard X-ray Burst Alert Telescope (BAT, 
used for the detection of GRBs, Barthelmy et al. (2005)), and the X-ray telescope (XRT, 
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used for their follow-up, Burrows et al. (2005)), are proving to be invaluable tools for sky 
monitoring and for observations of variable astrophysical objects such as novae. 



2.1. BAT Observations 

While awaiting GRBs, BAT normally operates in 'survey mode' in which spectral and 
imaging information is available with a time resolution which is typically 300s. A given part 
of the sky usually falls within the 1.4 sr BAT field of view many times per day. For example, 
64 observations of RS Oph, totalling more than 9 hours, are available in the period 3 days 
either side of the optical detection of the outburst. 

We have analysed all such data from t = —20 to +20 days using version 2.3 of the 
Swift software^. Standard filtering was used to remove data affected by bad quality, high 
background rates, source occultations, etc. Two non-standard programs were used for the 
overlaying of images and for the extraction of source count rates while allowing for the 
presence of strong sources in the field of view. The latter allows for the intensity of all 
sources in the field of view to be fitted simultaneously, along with detector background 
models. It was used to calculate weighted means of the estimated intensity at the position 
of RS Oph in 4 broad bands (14-25, 25-50, 50-100 and 100-190 keV). 

A clear detection of the source was made in the 14-25 keV band during the first 3 days 
from to (Figure 1). In an overlay of images formed from the corresponding data, standard 
software (batcelldetect) finds a 9.9o" source within 2.5 arc minutes of the position of RS Oph 
as the only unidentified source above 5o", confirming the reality of this detection. There are 
tantalising hints of emission in this same band in the preceding few days, and also a weak 
detection in the 25-50 keV band immediately following the outburst, but these should be 
treated with caution and their evaluation awaits a more detailed analysis. 

Standard software for spectral analysis was used. Spectra with 2 keV bins were extracted 
for each interval over which the pointing was unchanged (averaging 840 s). Systematic errors^ 
were taken into account and XSPEC^ Version 11.3.2 was used to fit models simultaneously 
to all the observations covering ~1 day. The signal to noise ratio was never good enough to 
justify models with multiple free parameters and so only the normalisation was fitted. 



^http://swift. gsfc.nasa.gov/docs/software/lhcasoft 
3per version 20051103 of the BAT CALDB 
''http://hcasarc/gsfc. nasa.gov/docs/xanadu/xspec 



-7- 



Unfortunately, due to a BAT reboot, no data are available at the time of the first 
XRT observation discussed below, and when data again became available the source was 
comparatively weak. We have fitted the BAT 14-50 keV spectra with the simple spectral 
models found using the XRT data as detailed below. Satisfactory x^ values were always found 
and allowing for the fact that BAT data are averages, and centered on slightly different times, 
the measured fluxes are consistent with extrapolations of the corresponding XRT spectra. 

The RXTE All-sky Monitor ^ observed the direction of RS Oph at several times during 
the period 2006 February 12-18 (t ~0 to 5 days). Little or nothing is detected in band a (1.3- 
3 keV) while in band h (3-5 keV), and more particularly in band c (5-12 keV), an increasing 
flux is seen. The fluxes are initially lower than expected from a model based on the flrst 
XRT observation, scaled to match the BAT data, but they reach the expected level by 
about February 15. This is consistent with an initially very high, but decreasing, absorption 
column. Sokoloski et al. (2006) have published an analysis of RXTE PCA observations 
starting from t = 3.1 days The hard X-ray light curve implied by their results is consistent 
with Figure 1, although the BAT observations extend higher in energy and the XRT data 
described below provide better low energy coverage and energy resolution below 10 keV. 



2.2. XRT Observations 

Table 1 gives details of the Swift XRT observations carried out over the flrst 26 days 
following the outburst. These began at t = 3.17 days (Bode et al. 2006; Bode 2006) and 
continued at a frequency which has given excellent temporal coverage. This means that 
the observations began much earlier in the evolution of the outburst, and we were able to 
follow it in much more temporal and spectral detail, than was the case for the EXOSAT 
observations in 1985. 

Swift software version 2.3 was used to process the data from the XRT. Source spectra 
were extracted from the cleaned Windowed Timing mode event-lists, using a box region of 
60x20 pixels (1 pixel = 2.36 arcsec); the same region was then offset from the source to obtain 
a background spectrum. During some of the observations the source was positioned over the 
bad CCD columns (Abbey et al. 2006, caused by a micrometeorite impact on 2005-05-27) 
and we corrected for the fractional loss of the PSF incurred when calculating the count rate 
and flux. 

Grade 0-2 events were chosen for the spectral analysis. The ftool xrtmkarf was used 
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to generate suitable ancillary response function files and these were used in conjunction with 
the most up-to-date response matrix (swxwt0to2_20010101v007.rmf). 

All spectra were grouped to a minimum of 20 counts per bin in order to facilitate the 
use of the x^ statistic in XSPEC, the results of which are discussed below. All fits included 
an energy scale offset that was allowed to vary between -0.1 and keV to account for 
possible CCD bias measurement uncertanties. Errors are given at the 90 per cent level (e.g., 
Ax^ = 2.7 for one degree of freedom). 

Figure 2 shows a sample of the resulting spectra at selected epochs and Table 1 gives 
further information, including the results of the model fits described below. The initial count 
rate at t = 3.17 days doubled in less than two days to peak at around 30 counts s~^, then 
fell again until the emergence of a new soft component beginning at Observation 8 (t = 26 
days) which dominated the XRT spectrum by Observation 9 (Osborne et al. 2006a, t = 29 
days) - see Figure 2. The gross evolution of the spectrum was a progressive softening during 
the period of these X-ray observations. 



3. Phases of Remnant Evolution 

In the earliest phase of the outburst, the ejecta from the WD surface are expected to 
traverse the binary system on a timescale t ~ 2d for Vej = 3000 km s~^ (Dobrzycka & Kenyon 
1994; Fekel et al. 2000). Subsequent to this, using models previously applied to the evolution 
of SNR, Bode & Kahn (1985) found that Phase I (where the ejecta were still important in 
supplying energy to the shocked stellar wind) ended in the first few days after the 1985 
outburst. Indeed, using the parameters derived by O'Brien et al. (1992) in Equations 5 and 
9 of Bode & Kahn (1985), we derive t ^ 6. 2d for the duration of Phase I (see Figures 3 and 
4). Thereafter, the remnant evolution can be well modelled by the instantaneous release of 
energy at a point. The Bode & Kahn (1985) analytical model of the 1985 outburst led to the 
conclusion that at the time the EXOSAT observations began (t = 55d), the remnant was in 
the transition between Phase II (where a blast wave is being driven into the wind and is so 
hot as to be effectively adiabatic) and Phase III (where the shocked material is well cooled 
by radiation). Initially, a double shock system is established, with the forward shock being 
driven into the stellar wind and a reverse shock being driven into the otherwise unshocked 
ejecta. In Phase II, the Primakoff similarity solution for this type of explosion into an r^^ 
density distribution (where r is the radial distance from the site of the explosion) gives 

r. = at2/3 fl) 
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where r^ is the shock radius and a is a parameter which is a function of the energy E 
communicated by the ejecta to the shocked wind, the mass loss rate in the red giant wind 
M and the red giant wind velocity u. Similarly, in Phase III, 

rs = be'\ (2) 

For a strong shock, the post-shock temperature Tg is given by 

where k is Boltzmann's constant and rfi = 10~^^ g is the mean particle mass, including 
electrons. 

The O'Brien et al. (1992) analytical models found good fits to the evolution of the rem- 
nant up to around t = 70 days after outburst when there was a significant decline in soft 
X-ray flux. This was explored in terms of numerical models of the forward shock running 
off the end of the red giant wind re-established between the 1967 and 1985 outbursts, as 
suggested by Mason et al. (1987). However, O'Brien et al. (1992) noted that discrepan- 
cies between the observed X-ray and models, particularly evident at low energies, could be 
explained by a contribution from unveiling of the central nuclear burning source (see be- 
low). The final EXOSAT observation on day 251 was proposed as being consistent with low 
luminosity remnant nuclear burning on the WD surface (Mason et al. 1987). 



4. Interpreting the Observations of the 2006 Outburst 

In the 2006 outburst, the Swift XBT observations reported here started towards the end 
of the period in which Bode & Kahn (1985) concluded that the reverse shock running into 
the ejecta may still be important (Phase I), then cover the period of the expected Phase I/II 
transition and extend into that of Phase II/III. 

From Table 1 and Figure 2 it can be seen that until around t = 26 days (Observation 8) 
the observed X-ray flux at first rose rapidly, then gradually declined, becoming progressively 
softer with time. However, as noted above, an excess at energies below around 0.7 keV was 
evident from t ~ 26d onwards and this varied dramatically thereafter (Osborne et al. 2006b). 
A fuller exploration of the evolution of the soft component will be given in Osborne et al. 
(2006c). 

Extensive model fitting of the early XRT results has been undertaken. The presence 
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of obvious line emission, both in the XRT spectrum and from a higher resolution Chandra 
X-ray observation on 2006 February 26 (Ness et al. 2006), strongly suggested that models 
of emission from a high temperature thermal plasma were most appropriate. Our first 
order approach has therefore been to fit a single temperature mekal model to the emission 
using XSPEC, yielding both line and continuum radiation, with the thermal bremsstrahlung 
contribution being increasingly dominant at higher temperatures. 

Up to and including Observation 7 at t = 18.17 days, the fit was performed over the 
whole 0.3 - 10 keV range of the XRT data. After this date, only data for E >0.7 keV were 
used in order to separate out emission from the additional, highly variable, soft component 
seen after this time. We have constrained the interstellar column to be [Nh]is = (2.4 ± 
0.6) X 10^^ cm~^ throughout (Hjellming et al. 1986) but we let the additional column from 
the overlying red giant wind, [A'"j:/]vi/, be a free parameter up to and including Observation 
7. [iV^/Jv^ was then fixed at 1.8 x 10^"*^ cm~^ thereafter as the new soft X-ray emission 
component began to dominate the emission at low energies. Spectral fits at several epochs 
using a solar abundance model are shown in Figure 2. We note that elemental abundance 
enhancements were required by Bode & Kahn (1985) in their model and also by Snijders 
(1987) and (Shore et al. 1996) from lUE observations of the 1985 outburst. However, we 
choose solar abundances until we have performed detailed analyses of the evolution of the 
shell (as described further in the Discussion section below). We also note that these single 
temperature fits, as may be reasonably expected, fail to reproduce all of the spectral detail 
and should not be overly interpreted at this stage. 

In the first three epochs, (Observations 1, 2a and 2b) the Fe K line was clearly detected 
at around 6.7 keV, as expected for the fitted plasma temperatures. In the first two epochs, 
there was also obvious excess emission at lower energies in the line compared to the single 
temperature mekal model fit. This could be characterised by an additional line at 6.4 keV 
which suggests X-ray refiection from a cold (or moderately ionised) gas as seen in some 
AGN and magnetic CVs (Beardmore et al. 1995; Done et al. 1995). We consider scattering 
in the red giant wind more likely in this case, a situation well documented in a number of 
supergiant X-ray binaries (e.g. van der Meer (2005) and references therein). The equivalent 
width of this line (~ 70 eV in Observation 2a) is consistent with fiuorescent emission by the 
surrounding gas with a column density ~ 6 x 10^^ cm~^ (Makashima 1985). 

Figures 3 and 4 illustrate the results for important parameters derived from the fits. For 
example. Figure 3 shows a monotonic decline in [iV^^jiy with time, in a manner qualitatively 
expected as the forward shock traverses the overlying wind. The forward shock velocity, Vg, 
has been derived from the best fit plasma temperature and Equation 3 above. It can also 
immediately be seen that the implied shock velocities are of the same order as those derived 
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spectroscopically (Bull 2006; Evans et al. 2006; Ness et al. 2006). The unabsorbed X-ray 
flux, Funahs) shows ail initial rise to a peak at around t = 5d (with Lx ~ 1.2 x 10^^ erg s~^ 
between 0.3 and 10 keV at d = 1.6 kpc) followed again by a subsequent decline. 

We can also undertake a more quantitative exploration of the basic model from these 
results in combination with VLBA observations obtained on 2006 February 26 (O'Brien et 
al. 2006a,b). These were taken almost simultaneously with Observation 5 of Swift (day 13.6), 
and showed a clumpy ring of emission which, if associated with the forward shock, would 
give a radius of approximately 2.1 x 10^^ cm sX d = 1.6 kpc. The value of Nh in the overlying 
wind (the net neutral column derived from the fits and shown in Table 1 and Figure 3) can 
be found from 

X MA 1 

where X is the mass fraction of hydrogen in the wind, itlh is the mass of a hydrogen 
atom, u is the wind speed and Tout is the outer wind radius (= wAt, where At is the time 
since the last outburst, ~ 21.04 years). With X = 0.7 for solar composition, -u = 20 km 
s~^ and M/u from O'Brien et al. (1992), this gives [iV/^]vi/ = 9 x 10^° cm~^, i.e. within a 
factor of two of the derived value from the single temperature model fits at this time (see 
Table 1). We note that [A'^ii^Jv^ derived in Equation 4 will be an underestimate for the value 
of the column for the whole shocked shell as it is calculated from the line-of-sight (front) 
edge. In addition. X-ray estimates of [A'^^jvi/ depend on the abundances in the wind (we 
have assumed these are solar) and furthermore Baumgartner & Mushotsky (2006) suggest 
that [Nh]is derived from X-ray observations may differ somewhat from 21 cm values as it 
also samples H2 in the ISM. 

After t ~ 6 days, (by which time the Bode & Kahn (1985) model expected Phase I to 
be effectively over) it can be seen from Figure 4 that Vs declines monotonically with time, 
following a power law decay from t ~ 10 days onwards. Here, Vg oc t~", where a ^ 0.6. This 
compares more precisely to the expected value of 0.5 for Phase III, rather than a = 0.33, 
expected in Phase II evolution (Bode & Kahn 1985). Equation 4 implies that [iVji/j^y oc t~" 
at early times (as Tout = 1.3 x 10^^ cm ^ r^), where a = 0.67 (Phase II) or 0.5 (Phase III) 
respectively. From Figure 3, it appears that after the implied end of Phase I, the data can 
in fact be better fitted by a = 0.5. Finally, after the expected end of Phase I, Funahs decayed 
as t~^'^, as expected from a well-cooled shock. 

The Sokoloski et al. (2006) observations of RS Oph with the RXTE PCA between days 
3 and 21 are also consistent with the basic shock model. Their results imply Phase II 
evolution after the first few days from the behaviour of their derived temperature, but are 
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more consistent with Phase III evolution during this period from the behaviour of the X-ray 
flux. The usable energy range of the PCA (2-25 keV) means however that their observations 
were less sensitive to N^ and they did not detect the emergence of the soft component. 

At t = 13.6 days, our fits to the Swift data imply Vg — 1700 km s^^. We note that we 
could obtain complete consistency between the VLBA imagery and simple derivation of Vg 
from the Swift model fits if a = 2.7 x 10^*^ cm s"^/^ and a larger distance oi d = 2.4 kpc, if the 
remnant were behaving precisely as one would expect for the Primakoff (Phase II) solution 
throughout. However, the time dependence of Vg, [Nh]w and Funabs are more like what one 
would expect in Phase III, but in this case h = 3.7 x 10^^ cm s~^/^ and for consistency with 
the VLBA results, the distance would have to be increased to 3 kpc. Such high values of 
the distance were accepted prior to 1985 (Bode 1987). On the other hand, Snijders (1987) 
notes that the failure to detect material at the typical velocities of the Carina Arm in lUE 
spectra places an upper limit on the distance to RS Oph of about 2 kpc. In addition of 
course, we have shown that there is a transition at around t = 6 days between remnant 
phases. Thus although the simple analysis is supportive of the qualitative correctness of the 
current models, as outlined below, detailed physical modelling is required to build a fully 
self-consistent picture. 



5. Conclusions and future work 

Our X-ray observations of the very earliest phases of the 2006 outburst of RS Oph are 
broadly consistent with the basic model of remnant evolution proposed by Bode & Kahn 
(1985) and further explored by O'Brien et al. (1992). In particular, it appears that Phase I 
may have ended by t ~ 6d as predicted. However, our first-order analysis suggests that the 
remnant moved rapidly into Phase HI in this outburst. 

Fits to the XRT data have also been attempted using vmekal where the abundances have 
been allowed to vary. For the single temperature fit, this resulted in a trend of decreasing 
elemental abundance enhancement with time, as would be expected as emission from the 
enriched nova ejecta at early times becomes increasingly less significant than that from the 
generally less enriched red giant wind. A series of fits using multi-temperature mekals was 
also performed producing improved fits to the data. However, detailed models involving 
full hydrodynamical simulations are ultimately required to fit the X-ray data in a physically 
meaningful way. These simulations produce an evolving range of temperatures and densities 
with radius from the central source (see O'Brien et al. (1992)). Resonant scattering of X-rays 
in the red giant wind may also be occurring, and therefore some of the X-ray emission line 
flux may be due to photoionization of this gas by the X-ray emitting shock. In addition. 
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VLBI results (O'Brien et al. 2006b) have shown that assumptions of spherical symmetry are 
called into question, and higher spatial dimension modelling may be required. We plan this 
more detailed modelling to be the subject of a future paper. 

In future papers we will report the results of our continued monitoring of the source 
with Swift. As well as attempting to understand the origin of the highly variable soft X- 
ray component (to be discussed in Osborne et al. (2006c), and which may possibly arise 
from continued nuclear burning near the Eddington luminosity on the WD surface), we 
are particularly interested in detailed investigation of three subsequent predicted phases of 
development: (a) the epoch when the forward shock is forecast to reach the end of the red 
giant wind at around t = 80 days, (b) the timescale for turn-off of the high luminosity phase 
of the TNR, and (c) the re-establishment of accretion in the central binary and the overlying 
wind from the giant. 

The authors are very grateful to the Swift Mission Operations Center staff for their 
superb efforts in supporting the observations reported here. We also acknowledge quick-look 
results provided by the ASM/RXTE team. P. Caraveo and A. de Luca contributed to work 
on the BAT data. SS gratefully acknowledges partial support from NSF and NASA grants 
to ASU. SPSE acknowledges the support of the UK's Nuffield Foundation. JPO, KLA, APB 
and MRG acknowledge support from PPARC J-UN is supported by the Chandra Fellowship 
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Table 1. Summary of RS Oph XRT Observations and Spectral Fits. 



Obs. 



Date (Day) 



Exposure 

(s) 



Count rate 
(count s^^) 



(keV) 



Correction 
Factor'^ 



[N^] 



w 



xVdoi 



'1022 cm-21 



Unabs. Flux'^ 
(erg cm~2 s~^) 



001 2006-02-16T00:05 (3.17) 

002a 2006-02-17T20:33 (5.03) 

002b 2006-02-21X00:10 (8.18) 

004 2006-02-23T19:37 (10.99) 

005 2006-02-26T10:18 (13.6) 

006 2006-02-28T10:31 (15.61) 

007 2006-03-03T00:04 (18.17) 
008^^ 2006-03-10T19:38 (25.99) 



664 


14.2 ± 0.2 


8.44l°|? 


1.2 


2.96 ± 0.15 


460/286 


999 


31.5 ± 0.2 


8.54+°i° 


1.1 


fi2+0-°3 

"•"^-0.02 


993/585 


844 


19.8 ± 0.2 


7 94+0.54 
'•^^-0.39 


1.0 


0.25 ± 0.02 


631/436 


950 


15.5 ± 0.2 


4.38 ± 0.23 


1.3 


0.25 ± 0.02 


484/327 


897 


12.3 ± 0.1 


3.30 ± 0.15 


1.1 


0.201H? 


465/289 


1041 


10.6 ± 0.1 


9 oo+O.lO 
^•O'-'-0.09 


1.0 


0.19 ± 0.02 


559/283 


738 


8.6 ± 0.1 


2 36+°-^^ 


1.0 


0.18_oo2 


353/206 


1605 


6.5 ±0.1 


-'-•"-'--0.02 


1.1 


0.18 


708/189 



2.0 X 10^9 

2.7 X 10-9 
1.3 X 10-9 
8.9 X 10-1° 

6.8 X 10-1° 

4.9 X 10-1° 
3.6 X 10-1° 

2.5x10-1° 



^Factor by which the count rate (and flux) has been corrected to account for bad columns on the CCD. This varies due to the 
changing placement of the source on the detector. 

I' Absorption above the ISM value of 2.4 x lO^i cm"^. 

■^0.7-10 keV. Bolometric corrections range from 1.5 to 1.7 for Observations 1 to 8. 

'iFits to 0.7 - 10 keV and additional column flxed at previous value (see text). 
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Fig. 1. — The flux from RS Opli in tlie bands 14-25 keV and 25-50 keV observed witli Swift 
BAT. Units are ergs cm"^ s~^ and are based on count rates in tliese bands assuming tlie 
spectral form found from the observation 001 XRT data. The circles show the flux in the 
same energy bands from extrapolation of the models in Table 1. Means have been calculated 
over 24 h periods and are represented by a point at the weighted mean time, with horizontal 
error bars stretching from the first to the last sample included in the mean . The very 
irregular sampling leads to uneven horizontal error bars. The dotted line is the time of the 
first optical detection of the outburst. 
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Fig. 2. — Selected Swift XRT spectra of RS Ophiuchi. Shown are data from t = days 3.17 
(black), 5.03 (red), 13.6 (blue), 25.99 (orange), and 29 (grey). The mekal, optically thin, 
collisionally-ionized plasma emission model fits shown are for days 3.17, 5.03 and 13.6, with 
their residuals plotted in the bottom panel. 
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Fig. 3. — Temporal behaviour of the overlying absorbing column in the red giant wind, 
[-^i?]vK) net of the assumed interstellar column of 2.4 x 10^^ cm~^, derived from the XRT 
spectral fits. A vertical dashed line on this Figure (and Figure 4) indicates the expected end 
of Phase I development. The diagonal dashed line is a power-law of index a = 0.5 . Note 
that Observation 8 is not shown as [iV//]vi/ was fixed at the value derived for Observation 7 
(see text for details). 
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Fig. 4. — The XRT-derived temporal behaviour of the shock velocity, Vs (points with error 
bars) and unabsorbed (0.7 - 10 keV) flux, Funabs (crosses). The diagonal dashed lines are 
power-laws with index a = 0.6 (for Vs points) and a = 1.5 (for Funabs - see text for details). 



